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The KOTO {K'^ at Tokai) experiment aims to observe the CP-violating rare decay 
Kl ^ TT^niy by using a long-lived neutral-kaon beam produced by the 30 GeV pro¬ 
ton beam at the Japan Proton Accelerator Research Complex. The Kl flux is an 
essential parameter for the measurement of the branching fraction. Three K^ neu¬ 
tral decay modes, Kl —>■ 37r°, Kl —>■ 27r°, and Kl —>■ 27 were used to measure the 
ATl flux in the beam line in the 2013 KOTO engineering run. A Monte Carlo simu¬ 
lation was used to estimate the detector acceptance for these decays. Agreement was 
found between the simulation model and the experimental data, and the remaining 
systematic uncertainty was estimated at the 1.4% level. The Kl flux was measured 
as (4.183 ± O.OlTstat. ±0.059sys.) x 10^ Kl per 2 x lO^"* protons on a 66-mm-long Au 
target. 
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1. Introduction 

The KOTO {K^ at Tokai) experiment is aimed at the first observation of the rare CP- 
violating decay of the long-lived neutral-kaon Kl —)• An important characteristic of 

this decay is its theoretical cleanness. The branching fraction (BF) is predicted to be 2.4 x 
10“^^ with a theoretical uncertainty of only 2.5% [T], in contrast to that of most other 
meson decays, especially in the B system, where the uncertainty can be as high as 10%. 
The study of —>• can also shed light on the origin of CP violation as it provides 

a direct precise measurement of the parameter rj |2], the height of the unitarity triangle, 
in the standard model (SM) of particle physics. In addition, —)• involves an s —)• d 

transition that is highly suppressed in the SM and thus is very effective in searching for new 
physics that modifies the SM flavor structure. The current experimental upper limit on the 
BF is 2.6 X 10“® at the 90% confidence level, obtained by the KEK-E391a experiment [3] at 
the KEK 12 GeV proton synchrotron (KEK-PS). 

The KOTO experiment [1] uses the experimental technique established by E391a. In order 
to improve the sensitivity by three orders of magnitude, a new dedicated beam line was 
constructed at the Japan Proton Accelerator Research Complex (J-PARC) [5j. The design 
value of the proton intensity delivered by the Main Ring (MR) at J-PARC is 100 times higher 
than that of the KEK-PS. To improve the background suppression, the main electromagnetic 
calorimeter and other detectors surrounding the decay volume were upgraded. 

The Kl flux generated by the proton beam is an essential parameter for the measurement 
of the BF of Kl —?■ The KOTO Collaboration measured the Kl flux and the momen¬ 

tum spectrum by reconstructing the charged decay modes iF/, —)• 7r+7r“7r° and Kl —)• vt+tt” 
in beam surveys with dedicated tracking detectors and calorimeters PE]. During data acqui¬ 
sition, the Kl flux is measured using the neutral decay modes Kl ^ Jtt®, Kl ^ 27r‘^, and 
Kl ^ 27, since the KOTO detector cannot measure the momentum of charged pions. In this 
paper, we report on KOTO’s first measurement of the Kl flux by using the three neutral 
decay modes with the main calorimeter, which was performed during an engineering run in 
January 2013. 

The paper is organized as follows. Section[^describes the experimental setup relevant to the 
measurement. Sectionillustrates the principle and the technique used in the measurement. 
Section reports on the Monte Carlo (MC) simulation used to measure the acceptance of 
each of the three modes. Section describes the analysis of the experimental data for each 
of the three modes. Section [^reports the results and their uncertainty estimations. Section]^ 
compares this to previous results. Section describes the conclusions. 

2. Experimental setup 

This section focuses on the experimental setup used for the measurement, including descrip¬ 
tions of the beam line and of the subsystems in the KOTO detector relevant for this 
analysis. 

2.1. Beam Line 

The KOTO detector is located at the Hadron Experimental Facility of J-PARC, at the end 
of a new 21-m-long neutral-kaon beam line (KL). A schematic top view of the KL beam 
line and the KOTO detector is shown in Fig. The protons are accelerated to 30 GeV by 
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Fig. 1: Schematic top view of the KL beam line used by the KOTO experiment. The primary 
proton beam comes in from the left. The KL beam line starts at the T1 target and the beam 
axis goes off at an angle of 16° with respect to the primary beam line. The front of the 
KOTO detector is located at the end of the beam line, about 21 m from the T1 target. 


the MR, extracted by using a slow extraction technique [8], and transported through the 
primary beam line to the facility [U] • The proton beam intensity is monitored by a secondary 
emission chamber (SEC), located after the extraction point. The primary proton beam, with 
a cross section of approximately 1 mm in radius, is injected into the production target (Tl). 
The target consists of a 66-mm-long gold target of 6 x 6 mm^ cross section. It is equally 
divided into six parts along its length with hve 0.2-mm-thick slits. 

The KL beam line is off-axis by an angle of 16° with respect to the primary proton 
beam line. The full kinematic reconstruction of the neutral pion in the Ki ^ decay 
requires a small diameter for the Kl beam and was achieved with the collimation scheme 
shown in Fig. The secondary particles produced at the target pass through the pair of 
collimators shown in Fig. to shape the beam. At the exit of the second collimator, the 
neutral beam has a square cross-section of 8.5 x 8.5 cm^ corresponding to a solid angle 
of 7.8 /rsr. Before entering the collimation region, the beam passes through a 7-cm-thick 
(12.5 Ao) lead absorber which removes most of the photons. A 2-T dipole magnet, located 
between the two collimators sweeps out charged particles. Short-lived particles decay in the 
long collimator region leaving only mesons, neutrons, and photons at the entry of the 
detector region. Table lists the composition and location of all the components along the 
KL beam line. The materials in the beam decrease the Kl flux by 60%, as estimated using 
a Geant4 mm based MC simulation. The collimation scheme of the KL beam line is 
described in Ref. [T^ . 

2.2. KOTO detector 

Figure shows a cross-sectional side view and the coordinate system of the KOTO detector. 
The signature of a —>■ decay is a pair of photons coming from the decay, without 
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Fig. 2: KOTO collimation scheme in the beam line coordinates |12j . 


Table 1; Composition and position of materials along the axis of the KL beam line: some 
components are part of the Kl.l charged-kaon beam line, which uses the same T1 target 
and is at an angle of 7° with respect to the KL beam line. The thickness is listed only for 
materials actually crossed by the beam. The regions between the Be vacuum window and the 
stainless steel vacuum window, and between the Kl.l front duct and the beam exit window 
are in a 2-Pa vacuum. The starting position is in beam line coordinates, whose origin is at 
the center of the T1 target. The KOTO detector starts at the front face of the Front Barrel; 
detector coordinates are measured with respect to this position. 


Name 

Material 

Thickness 

[mm] 

Starting position 
[mm] 

T1 target 

Au 

66 

- 

Vacuum window 

Be 

8 

247 

Vacuum window 

Stainless steel 

0.2 

3,097 

Photon absorber 

Pb 

70 

3,730 

Kl.l front duct 

Stainless steel 

0.2 

4,182 

Kl.l tail duct 

Stainless steel 

0.2 

5,510 

Collimator vacuum window 

Stainless steel 

0.1 

6,400 

1®* collimator 

Fe and W alloy 

- 

6,500 

2^^*^ collimator 

Fe and W alloy 

- 

15,000 

Beam exit vacuum window 

Polyimide 

0.125 

20,000 

Front Barrel 

- 

- 

21,507 

Csl calorimeter 


- 

27,655 


any other detectable particles. The energy and position of the two photons are measured with 
a cesium iodide electromagnetic calorimeter (Csl) [13]. Multiple charged-particle and photon 
detectors surround the decay volume to form a hermetic veto against any extra particles 
except neutrinos. The decay vertex of the Kl is reconstructed under the assumption that 
the two photons come from a vr® on the beam axis and that the vertices of the Kl and tt^ 
coincide. Finally the vr^ is required to have a large transverse momentum to balance the 
momentum carried by the two neutrinos. 
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Fig. 3: Cross-sectional side view of KOTO in detector coordinates, defined with respect to 
the front face of the Front Barrel (FB) detector. Ki particles come in from the left. The z- 
axis is along the beam direction, the y-axis points upward, and the x-axis goes into the page. 
The Csl calorimeter, located downstream of the decay volume, is surrounded by two photon 
detectors (FB and MB), charged-particle veto detectors (CV, LCV, BCV, and Hinemos), 
photon veto collar counters made of Csl crystals (NCC, CC03, CC04, CC05, and CC06) 
and forward charged-particle and photon detectors (BHCV and BHPV). The gap between 
the MB and the Csl calorimeter is filled with a photon detector (OEV). BHTS is a trigger 
scintillator which is used for detector performance checks. 


For this flux measurement, in addition to the Csl calorimeter, the Main Barrel 
(MB) [13] and the Charged Veto (CV) [15], which are described below, are used as veto 
detectors. 

The Csl calorimeter consists of 2716 undoped Csl crystals stacked in a cylindrical shape 
of 2-m diameter and 500-mm depth (27 Xq) along the beam direction. We use crystals of 
two sizes: 2240 small crystals of 25 x 25 mm^ cross section for the inner region and 476 
large crystals of 50 x 50 mm^ for the outer region. The calorimeter has a 148 x 148 mm^ 
square beam hole at its center. The Csl crystals have a typical light yield of 9 photoelec¬ 
trons (p.e.)/MeV, and an energy resolution (Te/E of 1.9%/-^/E^[Ge^ with a constant term 
of 0.6%. They are read out by photomultiplier tubes (PMTs) [T6] . 

The MB is a cylindrical photon veto detector 5.5 m long with 2-m inner diameter sur¬ 
rounding the decay volume. It consists of 32 modules, each made of 45 alternate layers of 
plastic scintillator and lead sheet. The inner 15 layers have 1-mm-thick lead sheets while 
the outer 30 layers have 2-mm-thick lead sheets. The thickness of the plastic scintillator is 
5 mm for both the inner and outer layers. The total thickness of the modules is 13.5 Vq. The 
scintillation light originating in the plastic scintillator is absorbed by wavelength-shifting 
fibers (WLSFs) and read out by PMTs that are connected at both ends of each inner and 
outer module. Figure]^ shows a cross-sectional end view of the Csl calorimeter and the MB 
photon veto detector. 

The CV is the main charged-particle veto detector. It consists of two layers, placed 5 cm 
and 30 cm upstream of the calorimeter. Each layer is made of 3-mm-thick and 69-mm-wide 
plastic scintillator strips of lengths between 490 and 1002 mm, assembled as shown in Eig.[^ 
The scintillation light is picked up by WLSFs and read out by multi pixel photon counters 
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Fig. 4: Cross-sectional end view of the Csl calorimeter and the MB photon veto detector. 


at both ends of each strip. The typical light yield of the strips is 186 p.e./MeV and the 
typical time resolution is 1.2 ns when averaging over the two ends m 
These three detectors are all contained inside a pressure vessel defining the decay volume, 
which is evacuated to 10“^ Pa to suppress interactions of the beam particles with the resid¬ 
ual gas. Due to a large amount of outgassing, the detectors inside the pressure vessel are 
separated from the high-vacuum region by a thin membrane and evacuated at a level of 
1 Pa. 

Analog signals from all detectors are digitized by ADC modules. A filter inside the ADC 
modules shapes the raw front-end signals into Gaussian pulses approximately 50 ns wide 
before digitizing them at a 125 MHz sampling rate and with 14 bit resolution [TH]. The 
digitized waveforms are transmitted to the data acquisition (DAQ) system via optical links 
and processed by a three-tiered trigger system m- The waveform consists of 64 consecutive 
8 -ns samples around the signal region, corresponding to a 512-ns time window. 

3. Kl flux measurement 

We performed an engineering run in January 2013 to check the performance of the KOTO 
detector, including the data readout and peripheral slow control systems. In addition to the 
commissioning of the detector, the beam flux with the detector in situ was measured. 
This section describes the principles of the Kl flux measurement, the run conditions under 
which the data were taken, and the triggers used during the run. 

3.1. Measurement principle 

We determined the Kl flux at the beam exit vacuum window (see Table. [^. The number of 
Kls were counted by using the neutral decay modes, Kl —)• Kl —)■ 27r^, and Kl —>• 27 , 
which have large statistics. Their branching fractions are listed in Table These decay 
modes have a common final-state configuration with all the decay products being photons. 
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Fig. 5: End view of the CV detector showing the length and orientation of the 69 mm-wide 
plastic scintillator strips in each of the two CV layers. 


Table 2: Branching fractions (BFs) of the three Kl neutral decay modes |20) used for the 
Kl flux measurement. 


Mode 

BF 

Kl SttO 

0.1952 ±0.0012 

Kl 2tt^ 

(8.64 ±0.06) X 10-4 

Kl 27 

(5.47 ±0.04) X 10-4 


They provide three independent results that can be compared and cross-checked against 
each other. 

The Kl ^ 2'y decays were reconstructed from events with only two photons assuming 
the nominal mass. The Kl —>• 27r'^ and Kl ^ decays were reconstructed from two 
orthogonal sets of events with exactly four and six photons detected in the Csl calorimeter, 
respectively. To reconstruct the vr^s, all possible combinations of two photon pairs were 
considered. The decay vertex position was calculated from the energy and position of the 
two photons in each pair, assuming the nominal vr® mass. The photon pairing was determined 
by minimizing the variance of the common vertex position, as described in Sect. |5.3.1 

The MB inner layers were used as a photon veto to suppress the —?• background 
in the —?■ 27r’^ and Ki ^ 2'y samples, when some of the photons fell outside the Csl 
fiducial area. The CV was used as a charged particle veto to remove backgrounds due to the 
Kl —)• 7r=*=e^z^ and —)■ 7r+7r“7r° decays by identifying charged-particles before they 

hit the Csl calorimeter. 
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Fig. Q: KI momentum distribution at the beam exit. Points with error bars are experimental 
data taken with a magnetic spectrometer [7]. The spectrum below 2000 MeV/c was measured 
with Kl —)• TT^TT”; that above 2000 MeV/c was measured with Kl —?■ vt+tt” and Kl —>• 
7r^7r“7r^. The red curve is a Gaussian fit of the data with a momentum-dependent resolution 
and it is used to generate the Kl spectrum in the MC simulation. The ht returns a peak 
momentum of 1.42 GeV/c. 

3.2. Run condition 

The data used for this analysis were accumulated during a 6-hour run taken on January 
16, 2013. The nominal beam power was 15 kW, corresponding to an intensity of 1.8 x 10^^ 
protons on target (POT) per spill. A spill had a duration of 2 s and a repetition rate of 6 s. 
The standard deviation of the beam intensity spill by spill was 0.7%. The total POT during 
the run, as measured by SEG, was (3.355 ± 0.013) x 10^®. The uncertainty was attributable 
mainly to the conversion factor from the SEC counts to POT. 

3.3. Trigger 

The Kl neutral decays used in the Kl flux measurement were triggered by requiring that 
the energy depositions in both the left and right halves of the Csl calorimeter be higher 
than a threshold. The common threshold for both halves of the calorimeter was 307.5 MeV, 
which resulted in a trigger rate of about 34 kHz. This trigger scheme is effective at selecting 
events with low missing energy. 

4. Monte Carlo simulation 

An MC simulation was used to estimate the acceptances of the detector for the Kl ^ 
37r°, Kl ^ 27r‘^, and Kl ^ 2'y decays. It was also used to evaluate the contamination of 
background events from other Kl decays. The Geant4 MC package (Geant4.9.5.p02) was 
used to simulate the geometry of the KOTO detector and the interactions between particles 
and materials inside the detector. Kl particles were generated at the beam exit with the 
momentum spectrum shown in Fig.[^ The spectrum was measured independently in a special 
run with a magnetic spectrometer [7j. 
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The detector response was corrected by using data taken in dedicated beam tests, bench 
tests with a ^^^Cs radioactive source, and cosmic-ray runs with the detector in situ. For the 
Csl calorimeter, we added to the MC simulation the measured light yield of each crystal 
and its nonuniformity along the depth of the crystal, the accuracy of the energy calibration, 
the pedestal fluctuations, the light propagation velocity along the depth of the crystal, and 
the time smearing due to the timing resolution. For the MB photon veto, the length of the 
detector resulted in a substantial dependence of the light yield and timing delay on the 
position of the hit along the module. This effect was simulated by using a light propagation 


velocity and light attenuation measured in data (see Sect. 5.4.2). For the CV, the light yield 


and timing resolution measured for each strip were used in the simulation. 


4 . 1 . Accidental overlay 

The overlap of accidental activity generated by simultaneous Kl decays or interactions of 
other beam particles was not simulated in the MC. Both of these accidental sources can result 
in additional energy deposition and possible distortion of the timing information recorded 
by the front-end electronics. To simulate such effects, we overlaid data collected with the 
accidental trigger described below to MC events passing the same selection criteria as the 
online trigger logic. The accidental data consisted of events triggered with the same trigger 


configuration described in Sect. 3.3, but recorded in a time window starting 512 ns or 768 ns 


before the actual trigger. The energy deposition measured in a randomly chosen accidental 
event was added to the energy deposition of a simulated event at the level of the ADC 
readout. If the energy of the overlapping event was larger that the energy of the simulated 
event, it determined the event timing at the detector. 


5. Analysis 

This section describes the calibration, the event reconstruction, and the event selection used 
for the Kl flux measurement. The photon, vr^, and Kl reconstruction methods are presented. 
A more detailed description can be found in Ref. |21j . 

5.1. Csl calibration 

The energy calibration of the Csl calorimeter was performed in three steps. An initial cal¬ 
ibration factor to convert ADC counts into energy was calculated by using vertical cosmic 
rays penetrating the calorimeter and assuming an energy deposit per unit path length in 
cesium iodide of 5.6 MeV/cm [20]. Next, from a large statistics Kl —)• 37r° sample, we tuned 
the calibration factor of each crystal by minimizing the overall width of the reconstructed 
Kl invariant mass distribution |22j . Finally, an absolute correction to the calibration factor 
was obtained by using events with of known decay positions. These data were collected 
in special runs with a 5 mm-thick aluminum plate placed at 3353 mm upstream of the Csl 
front face. Neutrons in the beam core interacted with the aluminium nuclei and generated 
7 r° which decayed immediately. The deposit energies of the two photons from the tt® —)• 27 
decays were measured by the Csl calorimeter. We scaled the absolute calibration factor of 
whole crystals to make the invariant mass of the two photons the same as the nominal 
mass. The invariant mass distribution of these events is shown in Fig. 
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Fig. 7: Reconstructed invariant mass of vr^ events collected in the aluminum plate runs. The 
distribution is htted using a Gaussian peak plus an exponential slope. The final fit is shown 
by the red curve. 



Fig. 8: Graphical view of the KOTO clustering algorithm: the colored boxes represent seed 
crystals. Gluster seeds with the same color share a cluster. The red circles show the cluster¬ 
growing radius of 71 mm. 


5.2. Clustering 

A photon hitting the Csl calorimeter induces an electromagnetic shower with a Moliere 
radius of 3.53 cm for cesium iodide [23], thus depositing energy over multiple crystals. First, 
we selected Csl crystals with energy deposits above 3 MeV and within 150 ns of each other. 
Such crystals were called “cluster seeds”. Next we grouped seeds into a single cluster if found 
within 71 mm of another seed, as shown in Fig. Once no other seed could be added, we 
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calculated the energy, coordinates, and timing of the cluster as follows: 


-^cluster — 



i=l 

^cluster — 

2^4=1 

V" P- ’ 

2^4=1 

2/cluster — 

V" P- ’ 

EILi 

^cluster — 



and 


(1) 

( 2 ) 

( 3 ) 

( 4 ) 


where n denotes the number of seeds in the cluster and e*, Xi, yi, and ti are the energy, 
position, and time of the i th seed, respectively. The time is normalized by an energy- 
dependent resolntion measured in test beam data as at- = -^/(57eip~+~(TM7^7eip~+~OT^ 
ns, with Cj in units of MeV m- 

Clusters with more than one crystal and an energy deposit over 20 MeV were called 
“photon clusters” and used in the vr® reconstruction. Single seed crystals that did not belong 
to any cluster were called “single-crystal hits” and used as a signature of accidentals in the 
Csl calorimeter. 

The center of energy (xduster, ^duster) was shifted from the actnal incident position of the 
photon dne to its nonzero incident angle on the surface of the Csl calorimeter, in conjunction 
with the depth of the main energy deposit. In addition, the energy of the cluster i?ciuster 
deviated from the true photon energy due to the shower leakage outside the calorimeter. 
We corrected for these effects by using correction maps extracted from an MC simulation in 
advance. The resulting position resolntion of the photon entrance position at the front face 
of the calorimeter was [GeV] mm with a constant term of 2.1 mm. 


5.3. Event reconstruction 

In the analysis of Kl —)• Kl —>■ 27r°, and ^ decays, the first step was to count 

the number of photon clusters. We required six or more photon clusters for Kl ^ Svr^, four 
or more photon clusters for Kl —)• 27r^, and two or more photon clusters for ^ 2'y events. 
Next, we selected the six, four, or two photon clusters closest in time and started the process 
of reconstructing the vertex and its four-momentum. 

5.3.1. Vertex reconstruction. For Kl ^ 2'y events, the event vertex was derived assuming 
that the two photon clusters have an invariant mass equal to the mass and originate 
from the beam axis at coordinates (0, 0, Zytx)- 
The reconstruction of the event vertex for Kl —>■ and Kl ^ 27r*^ events started by 

constraining pairs of photon clusters to have an invariant mass equal to the nominal 
mass. The calculation of the event vertex position (Wtx, Wtx, -^vtx) and timing Tvtx follows 
the same principles as what was done in the E391a experiment [3]; it is outlined in the 
following. 

Figure represents a graphical description of the reconstruction. Assuming that the 
invariant mass of the two photons is equal to the nominal mass [2^, the following 
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Fig. 9: Graphical description of the reconstruction: the two yellow lines show the paths 
of the photons from the event vertex to their entrance position at the front face of the 
calorimeter and 9 is their opening angle; are the energies deposited by the photons; 

2Csi is the z position of the front surface of the Csl calorimeter in detector coordinates; dz 
is the distance of the calorimeter from the event vertex; di, d 2 are the distances of the two 
photon hits from the event vertex; and ri, r 2 are their projection at the front face of the 
calorimeter. 


equation holds: 

= 2E-y^E^^{l — cos 9), (5) 

where E-y^ and E^^ are the energies of the photons and 9 is the angle between their directions. 
In addition, the following geometrical relations hold: 

ri 2 ^ = + d 2 ^ - 2 did 2 cosd, (6) 

di = + dz‘^, and (7) 

d2 = \fr^^^+~d^, ( 8 ) 

where ri 2 is the distance between the hit positions of the two photons at the front face of 
the calorimeter and dj (rj) is the distance between the hit position of the i th photon and 
the decay vertex (z-axis). From Eqs. , one can calculate the distance of the tt^ from 

the front face of the calorimeter dz, with an uncertainty derived from propagating the single 
photon energy and position resolutions. These in turn can be used to derive the position Z 
and uncertainty of the tt^ vertex in detector coordinates. 

The six (four) photon clusters in the Kl —)• [Kl —)• 27r‘^) samples can be paired to 
form three (two) in 15 (3) independent combinations. To find the correct photon pairing, 


we minimized the “pairing variance” Xz- 


2 ^ 

W 2. ’ 

i=l 

(9) 

ELi 

EEi ’ 

(10) 


where n is the number of pairs in a given combination. We assigned the Z of the combination 
to the smallest Xz as the Kl decay position Zvtx- 
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The Kl vertex position in the x-y plane was calculated as the interpolated position between 
the T1 target and the center of energy on the Csl front face at the position assuming 
that the T1 target is a point source. 

Once the three-dimensional vertex position was determined, the generated time of each 
photon at the vertex Tj was calculated from the photon cluster time tj after correcting for 
the time-of-flight as Tj = tj — dj /c, where dj is the distance between the event vertex and 
the calorimeter front face (see Fig. [^, and c is the speed of light. Finally, the time of the 
Kl decay time Tvtx was defined as the weighted mean of the photon generated times: 


Tvtx — 


ET=iT,/al 
ESi i/^l ’ 


( 11 ) 


where at - is the timing resolution of a cluster given by the sum in quadrature of the energy- 
dependent term, 3.8/^/£|)~[MeVj ns, and the constant term 0.19 ns. 


5.3.2. Mass and momentum reconstruction. So far we have assumed that the photon 
pairs in Kl ^ and Kl ^ 27r^ decays have the nominal mass. With the reconstructed 
vertex position in hand, these constraints were removed and the four-momenta of the tt® were 
recalculated by assuming the vertex position (Xvtxj Wtx, -^vtx)- The four-momentum of the 
initial Kr was obtained by summing over the four-momenta of the vr*^. For the Kl ^ 2"f 
decay, the four-momentum of was calculated from the four-momenta of the two photons 
assuming the vertex position ( 0 , 0 , .^vtx)- 


5.4. Event selection 

The selection of —?• events by requiring exactly six photon clusters in the Csl 

calorimeter resulted in a sample with 10% background contamination in the MC simulation. 
A similar requirement of exactly four photon clusters for Ki ^ 2 tt^ events, and exactly two 
for iFi —)• 27 events, resulted in signal-to-background ratios of 0.18% and 0.25%, respec¬ 
tively. Selection criteria (cuts) on the kinematics of the event and in the presence of extra 
particles were applied for each mode to improve the signal-to-background ratio, as described 
in Sects. [T4.1 and 5.4.2 Most cuts are common to all modes, as summarized in Tables Hi 
andlH 


5.4.1. Kinematic cuts, 
o ATvtx cut 

The difference between the time Tj of each photon in Eq. ( |11[ ) and the reconstructed 
vertex time T^tx was required to be less than 3 ns. This cut reduced the contamination 
of accidentals in signal events. 

® (A^miii) kmin) cut 

The position of the innermost photon was required to be outside an area of 120 mm 
X 120 mm from the Csl calorimeter center. This cut removed photons whose shower 
leaked into the beam hole, 
o .Rmax cut 

The position of the outermost photon was required to be inside a radius of 850 mm 
from the center of the Csl calorimeter. This cut removed photons whose shower leaked 
outside the calorimeter fiducial volume. In the ATl —>■ 27 analysis, the outermost photon 
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was also required to have a minimum radius of 450 mm in order to remove background 
events from the ^ mode decaying near the front face of the Csl calorimeter. 

O Erain CUt 

The energy of each photon was required to be larger than 50 MeV. This cut removed 
photon clusters with poor energy and position resolutions. 

O drain CUt 

The distance between the hit position of any two photons at the front face of the Csl 
calorimeter was required to be larger than 150 mm. This cut reduces the probability 
of misreconstructing a single photon into multiple clusters or multiple photons into a 
single cluster. 

° xLape cut 

In order to distinguish electromagnetic showers generated by single photons from showers 
generated by multiple photons or hadronic interactions, we defined the cluster shape 
variable: 


Xshape 



( 12 ) 


where N is the number of crystals involved in the cluster, denotes the observed 

(expected) energy deposit in each crystal, and is the uncertainty on the expected 
energy deposit. The single-crystal energy deposit, and its uncertainty, were derived from 
an MC simulation by using photons of different incident energies and polar and azimuthal 
angles [7j. A cut of Xshape ^ ^ used for the ATl —)• 27r° and —)• 27 events, 
o AM^o cut 

The reconstructed invariant mass of the was required to be within 10 MeV/c^ of the 
nominal for the —)• mode and within 6 MeV/c^ for the AT/, —?■ 27r*^ mode. 

This cut rejected background events with mispaired photon clusters, 
o AMkl cut 

The reconstructed invariant mass for Kl —)• 2tt^ events was required to be within 
15 MeV/c^ of the nominal Kl mass (20]. This cut was effective for removing Stt^ back¬ 
ground events with two photons not reconstructed inside the Csl calorimeter fiducial 
volume, 
o Pt{Kl) cut 

The reconstructed Kl transverse momentum for the Kl —?• 27 decay mode was required 
to be less than 50 MeV/c, higher than the 30-MeV/c peak observed in the Pt{Kl) 
distribution of simulated events. This cut rejected background events, in particular from 
the Kl —?■ mode, where some missing transverse momentum was carried away by 

undetected particles, 
o Zvtx cut 

The reconstructed Kl vertex was bounded to be between 2000 mm and 5400 mm in the 
detector coordinate for the Kl —>• and the Kl —?• 27r’^ modes, and between 2000 mm 

and 5000 mm for the AT^ —?■ 27 mode. 

^ (A^coe? h^oe) cut 

The event center-of-energy, given by the energy-weighted x and y coordinates of all 
photons at the front face of the calorimeter, was required to be within a 60 x 60 mm^ 
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area around the center of the calorimeter. This cut helped to reject background events 
with missing energy coming from undetected particles in the Kl ^ 27r*^ mode, 
o T,Eii2 cut 

The sum of the energy deposited by photons in both the left and right halves of the Csl 
calorimeter was required to be above 350 MeV, above the trigger threshold of 307.5 MeV. 

5.4-2. Veto cuts. 
o Csl veto 

In addition to measuring energies and positions of photons, the Csl calorimeter also 
served as a veto detector for extra photons. Events were removed if one or more 
photon clusters were present in addition to the six, four, and two photons for the 
Kl ^ 37r‘^, Kl —>• 27r‘^, and Kl ^ 2"/ modes, respectively. The additional photon clus¬ 
ter was counted only if its timing was within ±10 ns of the average time of the photon 
clusters used in the event reconstruction, 
o Csl single-crystal hit veto 

Events with “single-crystal hits” not included in any cluster were rejected. Eluctuations 
in the electromagnetic shower occasionally result in single crystal hits near a photon 
cluster. To avoid signal acceptance loss, the minimum energy threshold for a single¬ 
crystal hit was determined as a function of the distance from the closest photon cluster, 
d, as: 


Ethr = (23.4 — 0.034 d[mm]) MeV 
Ethr = 3 MeV 


tor iUU < a < DUU mm. 


for d > 600 mm. 


Vxo; 


(14) 


Furthermore, the time of the single-crystal hit was required to be within ±10 ns of the 
closest photon cluster. We did not veto events with any activity within 100 mm of the 
photon clusters. 

CV veto 

Penetrating charged particles deposit around 0.5 MeV in the CV. Based on this condition 


and the typical CV timing resolution (see Sect. 2.2), we defined a “CV hit” as any activity 


in a strip with a minimum energy deposit of 0.4 MeV in a time window of ±10 ns from 
the event vertex Tvtxj after correcting for particle time-of-flight. The energy and timing 
of the hit were calculated from the sum of the energy deposit and the average of the time 
measured at the two ends of the CV strip, respectively. By rejecting events with any 
hits in the CV, background events with charged particles were effectively suppressed. 
Inner MB veto 

The role of the inner MB veto was to veto the Kl —)• 37r^ background for the Kl —>■ 27r*^ 
and Kl ^ 2'y modes. The energy threshold was set at 5 MeV, and the time window, 
corrected for time-of-flight analogously to what was done for the CV, covered the interval 
[—26, ±34] ns. The inner MB timing resolution was 1.3 ns for a 5-MeV energy deposit. 
A wide time window was chosen to allow for time differences between the MB and the 
Csl calorimeter due to the MB length (5.5 m) along the beam direction. The energy, 
position, and time of a hit within the MB, after correcting for attenuation effects and 
the propagation delay in the WLSFs, were calculated using the signals from both the 
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Table 3: Cuts used in the selection for the —)• mode 


Kinematic cut 

Min. 

Max. 

ATytx 


3 ns 

(^min? ^^in) 

120 mm 


^max 


850 mm 

-^min 

50 MeV 


^min 

150 mm 


AM^o 


10 MeV/c^ 

AMk, 


15 MeV/c^ 

Zvtx 

2000 mm 

5400 mm 


350 MeV 


Veto cut 

Energy thresh. 

Time window 

Csl 

20 MeV 

— 10 ns —hlO ns 


upstream and downstream ends of a MB module as: 

^up T ^down 


^MB = 
ZUB = 
-£^MB = 


^up ^down 
2 


^ '^prop; and 


-up 


+ 


6down 


exp [-zmb/ (A + azMs)] exp [zmb/ (A - o^mb)] ’ 


(15) 

(16) 
(17) 


where eup(edown) is the visible energy and tup(tdown) the hit time at the upstream (down¬ 
stream) side, Zmb is the hit position along the beam axis with respect to the center of 
the MB module, Uprop = 168.1 mm/ns is the propagation velocity of light in the WLSF, 
and A = 4923 mm and a = 0.495 parameterize the effects of the attenuation of the signal 
along the length of the MB. 


6. Results 

In this section we derive the result of the flux from the measured Ki, yield in each of 
the Kj, —7- 37r°, —?■ 27r®, and Ki, —)• samples, after normalizing to the known number 

of protons on target. The Kj, yield V is defined as T = A^^'^^*^/(AxBF) where is the 

number of reconstructed data events after applying all the selection cuts and the background 
subtraction that is described in this section. A is the product of the decay probability and 
the acceptance calculated as the ratio of MC events passing the selection cuts over the total 
number of Kj, that passed the beam exit and decayed to the given decay mode, and BF is 
the branching fraction for the specific decay mode listed in TableIn the following sections, 
the inputs to the yield measurement for each of the three modes are presented. 


6.1. yield for Ki —)■ Svr^ mode 

The reconstructed invariant mass distribution for Ki ^ 3'7r^ candidate events after applying 
all the selection cuts in Table except for the cut, is shown in Fig. 10, The MC 

distributions in the figure are normalized so that the number of reconstructed events after 
applying all the selection cuts in the MC simulation with the background contamination is 
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Table 4: Cuts used in 

the selection for the KTl —?■ 2 Tr^ mode 

Kinematic cut 

Min. 

Max. 

ATvtx 


3 ns 

(-^min? ^^in) 

120 mm 


-Rmax 


850 mm 


50 MeV 


^min 

150 mm 


^shape 


5 

AM^o 


6 MeV/c^ 

AMk, 


15 MeV/c^ 

■2^vtx 

2000 mm 

5400 mm 

(-^coe? ^oe) 


60 mm 


350 MeV 


Veto cut 

Energy thresh. 

Time window 

Csl 

20 MeV 

— 10 ns —hlO ns 

Csl single-crystal hit 

3-20 MeV 

— 10 ns —hlO ns 

CV 

0.4 MeV 

— 10 ns —hlO ns 

Inner MB 

5 MeV 

—26 ns —1-34 ns 


Table 5: Cuts used in 

the selection for the Kl ^ 2 'y mode 

Kinematic cut 

Min. 

Max. 

ATvtx 


3 ns 

(-^min? ^^in) 

120 mm 


-Rmax 

450 mm 

850 mm 

^shape 


5 

Pt{Kl) 


50 MeV/c 

•^vtx 

2000 mm 

5000 mm 


350 MeV 


Veto cut 

Energy thresh. 

Time window 

Csl 

20 MeV 

— 10 ns —hlO ns 

Csl single-crystal hit 

3-20 MeV 

— 10 ns —hlO ns 

CV 

0.4 MeV 

— 10 ns —hlO ns 

Inner MB 

5 MeV 

—26 ns —1-34 ns 


equal to that in the data. The background contamination was estimated to be less than 0.1%, 
smaller than both the statistical uncertainty (0.39%) and the uncertainty in the nominal BF 
(0.61% [2Q]). The final number of reconstructed events in the data and the MC acceptance 
are listed in Table The Ki yield at the beam exit for the Kl —?■ mode was estimated 

to be (7.835 ± 0.033) x 10®. The uncertainty includes only the statistical uncertainties in the 
data and MC inputs. 


18 




















400 420 440 460 480 500 520 540 560 580 600 


Reconstructed Kl Mass [MeV/c^] 

Fig. 10: Reconstructed Kl mass spectrum in the Kl —>• data sample. The black points 

represent the experimental data, the open red histogram represents the Kl ^ Stt^ MC sim¬ 
ulation, and the solid blue histogram represents the estimated contribution from the other 
Kl decay modes. The lower plot shows the ratio between the experimental data and the 
MC simulation. The error bars reflect only the statistical uncertainties. 


Table 6: Number of reconstructed events in the data after applying all the selection cuts for 
the Kl —)• Svr® mode and subtracting the background contamination. The MC acceptance 
for the signal and the dominant background contribution, before and after correcting for 
their BFs are also listed. The uncertainty in the acceptance reflects only the MC statistics. 
The last column shows the relative contribution of MC events for each decay mode. 


Ardata 
^ '^rec 

65, 284 ± 256 




A 

HxBF 

(%) 

Kl ^ 37r0 

(4.269 ±0.007) X 10"“^ 

(8.332 ±0.013) X 10-5 

99.95 

+ 

1 

o 

(1.74 ±0.19) X 10"^ 

(2.18 ±0.23) X 10"® 

0.03 


6.2. Kl yield for Kl ^ 2%^ mode 

The reconstructed invariant mass distribution for Kl ^ 27r^ candidate events after applying 
all the selection cuts in Table except for the AMk^ cut, is shown in Fig. 11, The MC 
distributions in the figure are normalized to the yield measured for the Kl —)• mode. 
The background contamination was estimated by minimizing: 

N 


x' = Yl 


(di - asi - (dbiY 


^ (o-d,y + 


(18) 


with respect to the signal and background scale factors, a and /3. Here, N is the number of 
bins, di ± ad- is the number of data events in the i th bin and its statistical uncertainty, and 
Si ± asi {bi ± (TbJ is the number of MC signal (background) events and their relative uncer¬ 
tainties normalized to the Kl —?■ flux result. The minimization returned a = 0.993 ± 

0.034 and /3 = 0.956 ± 0.036 with x^/ndf (number of degrees of freedom) = 53.3/48. The 
calculated number of background events was 120.8 ± 6.2 out of 1129 total Kl 
candidate events. 
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Fig. 11: Reconstructed Kl mass spectrum in the Kl —)• 27r° selected data sample. The black 
points represent the experimental data, the open red histogram represents the Kl —>■ 27r*^ 
MC simulation, the light-blue histogram represents the estimated contribution from the 
Kl ^ Stt^ mode, and the solid blue histogram represents the contribution from the other 
Kl decay modes. The lower plot shows the ratio between the experimental data and the 
MC simulation. The error bars reflect only the statistical uncertainties. 


Table 7: Number of reconstructed events in the data after applying all the selection cuts for 
the Kl ^ ‘i-K^ mode and subtracting the background contamination. The MC acceptances 
for the signal and the dominant background contribution, before and after correcting for 
their respective BFs are also listed. The uncertainty in the acceptance reflects only the 
MC statistics. Only the backgrounds contributing more than 0.1% are listed; their relative 
contributions are shown in the last column. 


Tydata 

'^rec 

1008.2 ± 34.2 




A 

AxBF 

(%) 

Kl 27r° 

(1.487 ±0.012) X 10"® 

(1.284 ±0.011) X 10"® 

88.8 

Kl 37r0 

(8.23 ±0.29) X 10"^ 

(1.606 ±0.056) X 10"'^ 

11.1 


The number of reconstructed signal events and the MC acceptance are summarized in 
Table The Kl yield at the beam exit using the Kl ^ 2tt^ mode was estimated to be 
(7.85 ± 0.27) X 10®. The uncertainty includes only the statistical uncertainties in the data 
and MC inputs. 


6.3. Kl yield for Kl ^ 2^ mode 

The reconstructed transverse-momentum distribution for Kl ^ 2y candidate events after 
applying all the selection cuts in Table[^ except for the Pt{Kl) cut, is shown in Fig. 12 The 
MC distributions in the figure are normalized to the yield measured for the Kl —>• mode. 

Analogously to what was done in the Kl ^ 27r^ analysis, the background contamination in 
the candidate event sample was extracted by minimizing the in E< 1 - (18). The minimization 
returned a = 0.973 ± 0.025 and (3 = 1.023 ± 0.031 with x^/ndf = 67.1/78. The calculated 
number of background events was 200.9 ± 8.4 out of 1890 total Kl —>■ 27 candidate events. 
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Fig. 12: Reconstructed Kl transverse-momentum spectrum in the Kl —>■ 27 selected data 
sample. The black points represent the experimental data, the open red histogram represents 
the iFi —)• 27 MC simulation, and the light-blue, green, and orange histograms represent the 
estimated contribution from the Kl ^ 27r°, Kl ^ and iFgS decay modes, respectively. 
The blue histogram represents the other Kl decay modes. The lower plot shows the ratio 
between the experimental data and the MC simulation. The error bars reflect only the 
statistical uncertainties. 


Table 8 : Number of reconstructed events in the data after applying all the selection cuts for 
the Kl —)• 27 mode and subtracting the background contamination. The MC acceptances for 
the signal and the dominant background contributions, before and after correcting for their 
respective BFs are also listed. The uncertainty in the acceptance reflects only the MC statis¬ 
tics. Only backgrounds contributing more than 0.1% are listed; their relative contributions 
are shown in the last column. 


^data 

''rec 

1689.1 ± 44.3 




A 

TxBF 

(%) 

Kl 27 

(4.035 ± 0.020) X 10"® 

(2.207 ±0.011) X 10"® 

89.8 

Kl 37r° 

(1.133 ±0.034) X 10-® 

(2.212 ± 0.066) X 10"'^ 

9.0 

Kl 27r° 

(2.58 ±0.16) X 10"® 

(2.23 ±0.14) X 10"® 

0.9 

Kl —>■ vr^e^z^ 

(1.4 ±0.4) X 10"® 

(5.68 ±2.15) X 10"^ 

0.2 


The number of reconstructed signal events and the MC acceptances are summarized in 
Table The Kl yield at the beam exit using the Kl —)• 27 mode was estimated as (7.65 ± 
0.20) X 10®. The uncertainty includes only the statistical uncertainties in the data and MC 
inputs. 

6.4- Kl flux 

The measured Kl yield for each decay mode is summarized in Table The three results are 
consistent within the statistical uncertainties. Their weighted average gives the final Kl flux 
of l^ve = (7.831 ± 0.032) X 10®. The uncertainty includes only the statistical uncertainty. 
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Table 9: Kl yield obtained for the Kl —)• Stt*’, Kl —)• 27r°, Kl —)• 27 decays, and their 
weighted average. The last column gives the relative yield for each decay mode with respect 
to the average. 


Mode 


Kl yield (xlO®) 

Relative yield 

Kl^ 

0 

CO 

7.835 ± 0.033 

1.001 ± 0.004 

Kl^ 

to 

0 

7.85 ±0.27 

1.002 ± 0.034 

Kl^ 

27 

7.65 ±0.20 

0.977 ± 0.026 

Average 

7.831 ± 0.032 



The absolute Kl yield in units of 2 x 10^“^ POT, which is the J-PARC design value for the 
number of protons on target per spill, was calculated from the Kl yield as: 

Number of Kl _ Tave (2 x 10^^) 

2 X 10^^ POT etrg A^pot 

= (4.188 ± 0.017) X 10^ (19) 

where etrg = 11.15% is the KOTO trigger efficiency corrected for losses due to the DAQ 
dead time, and A'poT is the total number of protons on target collected during the run. The 
uncertainty includes only the statistical uncertainty of the three modes. 


6.5. Systematic uncertainties 

In this subsection, we describe the sources and estimation methods for the systematic 
uncertainties that affect the ffux measurement. 

The uncertainty in the signal acceptance coming from the modeling of the calorimeter was 
estimated by comparing the effectiveness of each kinematic and Csl veto cut in data and 
MC. We calculated the single-cut fractional difference, defined as the ratio of data and MC 
efficiency for a given cut after all others have been applied. Figure 13 shows the single-cut 
efficiency of data and MC and their fractional differences for all the Csl-based cuts in the 
Kl ^ analysis, including the extra cluster veto cut. By summing in quadrature all the 
fractional differences, we obtained a total systematic uncertainty coming from the modeling 
of the Csl calorimeter of 1.38%. The systematic uncertainties for the other two modes were 
calculated in an analogous way and found to be 2.18% and 3.90% for the Kl —>• 27r° analysis 
and the Kl ^ 2^ analysis, respectively. Figures 14 and 15 show the single-cut efficiencies 
and their fractional differences for Kl ^ 27r^ and Kl —?■ 27, respectively. 

The systematic uncertainty in the acceptance from the modeling of the veto detectors’ 
response has two components: uncertainty from accidental losses and from backsplash losses. 
The accidental loss was due to accidental activity depositing energy in any veto detector 
in coincidence with a real photon in the Csl calorimeter. The backsplash loss was caused 
by particles belonging to a photon electromagnetic shower escaping the front of the Csl 
calorimeter and generating secondary activity in the veto detectors. The uncertainty in the 
modeling of both sources was studied by comparing the efficiency of the veto detector cuts 
in data and MC for Kl ^ 3'7r^ events, for which no CV or inner MB veto cuts were applied. 
The change in yield was measured to be 1.00% after applying the inner MB veto cut and 
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Fig. 13: Single-cut efficiency (top) and the fractional difference (bottom) for each Csl 
selection cut, including the extra cluster veto cut, used in the Kl ^ analysis. 



Fig. 14: Single-cut efficiency (top) and the fractional difference (bottom) for each Csl 
selection cut, including the extra cluster veto cut, used in the Kl —^ analysis. 


0.65% after applying the CV cut. These differences were added as systematic uncertainties 
in the Kl ^ 21 ^^ and Kl ^ 2 ’j selections. 

Other sources of systematic uncertainties were estimated by changing the offline energy 
threshold from 350 MeV to the online trigger threshold value of 307.5 MeV. The resulting 
change in the number of selected events was taken as the systematic uncertainty of the SFii/2 
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Fig. 15: Single-cut efficiency (top) and the fractional difference (bottom) for each Csl 
selection cut, including the extra cluster veto cut, used in the Kl —>• analysis. 


cut. The Kl yield calculation used the Particle Data Group (PDG) branching fraction central 
values for the three neutral decay modes; the uncertainties on the central values reported in 
Table have been considered as a source of systematic uncertainty. Finally, the conversion 


factor from counts in the SEC to proton intensity used in Sect. 3.2 originates the 0.39% 
uncertainty in the A^poT which is common to all the decay modes. 


All the sources of systematic uncertainties are summarized in Table 10 , The largest source 
is the modeling of the Csl calorimeter. All the other uncertainties are smaller than the 
statistical uncertainty for a given mode. The final uncertainty has been calculated by adding 
in quadrature all the statistical and mode-dependent systematic uncertainties for a single 
mode, taking their weighted average, and adding in quadrature the mode-independent A^poT 
uncertainty. From Eq. (19), the final Kl flux result was (4.183 ± 0.017stat. i 0.059sys.) x 10^ 


Kl per 2 x 10^^ protons on target. Figure [I^ compares the Kl flux result separately for the 
three modes. 


7. Discussion 

The result of this paper can be compared to the previous Kl flux measurement of (4.19 ± 
0.09](q 44) X 10^ per 2 x 10^^ POT, obtained from data taken in a dedicated beam survey 
run in February 2010 [6]. The experimental running conditions during the two measurements 


are summarized in Table 11 Although the T1 targets were made of different materials (An 


in 2013 vs. Pt in 2010), they have the same proton interaction length Xp within 3% 

Scaling the Kl flux measured here to the J-PARC design values of 300 kW for the beam 
power and 0.7-s spill duration every 3.3 s [Ij, we predict a, Kl flux 100 times larger than 
that available in the previous Kl —>• search experiment. Together with the upgrades to 

the detector, the experiment should reach the sensitivity of the standard model prediction 
for the Kl — t 'k^vV search over a period of three Snowmass years (3 x 10^ s). 
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Table 10: Summary of the systematic uncertainties. The first five sources are added in 
quadrature to obtain the mode-dependent systematic uncertainty of each decay channel. 
The A^pot uncertainty is independent of the three modes. 


Source 

Kl 37r0 

Kl 27r0 

Kl 27 

Csl calorimeter modeling 

1.38% 

2.18% 

3.90% 

Main Barrel modeling 

- 

1.00% 

1.00% 

Charged Veto modeling 

- 

0.65% 

0.65% 

SKi /2 cut 

0.24% 

0.41% 

0.04% 

PDG branching fraction jjnj 

0.61% 

0.69% 

0.73% 

Mode-dependent 

1.53% 

2.61% 

4.14% 

Mode-independent (Vpqt) 


0.39% 



XlO^ 



Fig. 16: Kl flux result for each of the Kl —>■ Kl —)• 27r°, and Kl —)• 27 modes. The 

hatched area shows the statistical uncertainty and the vertical bar shows the quadratic sum 
of the statistical and mode-dependent systematic uncertainties. The horizontal dashed line 
shows the weighted average value. 


8. Conclusion 

In this paper, we have described the Kl flux measurement for the KOTO experiment with 
data taken during a detector commissioning run in January 2013. The measurement was 
done by using three Kl neutral decay modes: Kl —>■ Kl 2tt^, and Kl —>■ 27. The 

results for the three decay modes agreed with each other within the statistical uncertainties. 
Systematic uncertainties were estimated based on the reproducibility of the data selection 
efficiency in the Monte Carlo simulation. The final Kl flux was (4.183 ± 0.017 ± 0.059) x 10^ 
per 2 X 10^^ protons on target, where the first uncertainty was statistical and the second was 
systematic. This result is in agreement with a previous measurement done by the KOTO 
Collaboration during a dedicated beam survey run in February 2010. 
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Table 11; Comparison of target material and target thickness, both in units of mm and 
proton interaction length Xp, for the two KOTO Kl flux measurements. The last two 
columns summarize the typical beam power conditions and the Kl decay modes used for 
each measurement. 


Period 

material 

Target 

thickness 

Xp 

MR beam power 

Measured mode 

Feb. 2010 

Pt 

60.0 mm 

0.658 

1 kW, 1.5 kW 

Kp — )■ 7r^7r“7r'^ 






O 

CO 

t 

Jan. 2013 

Au 

66.0 mm 

0.640 

15 kW 

Kl 27r0 


(1 

-mm slits included) 


Kl 27 
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